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NOT INTEGRABLE 
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Abstract 

The perturbation of an exact solution exhibits a movable transcendental essential 
singularity, thus proving the nonintegrability. Then, all possible exact particular 
solutions which may be written in closed form are isolated with the perturbative 
Painleve test; this proves the inexistence of any vacuum solution other than the 
three known ones. 
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1 Introduction 



The Bianchi IX cosmological model ||, |33|, |, || |, |, [L0|, 0, |T| is governed by 
a system of three coupled second order differential equations 

a 2 (Log A)" = A 2 -{B- Cf and cyclically, (I) 

where a prime means a derivation with respect to the "logarithmic time" r, and a 2 
is 1 or —1 according as the metric is asymptotically Euclidean or Minskowskian. An 
equivalent definition is 

cr 2 (Logc<Ji)" = cul + tul — uj 2 uj 2 /uj 2 and cyclically, (2) 

with 

A = UJ2UJ3 1 oj\ , uj\ = and cyclically. (3) 
This system possesses the first integral 

/ = a 2 {(LogBy(LogCy + (LogCy(LogAy + (LogAy(LogB)'} 

+A 2 + B 2 + C 2 - 2(BC + CA + AB) (4) 
= 2o- 2 {(Log^ 2 )'(Logu;3) / + (Logu^^Logu;!)' + (Logw^^Log^)'} 
^{(Logc^)' 2 + (Log^ 2 )' 2 + (Log^ 3 )' 2 } 

,-2 2 2i -2 2 2 1 -2 2 2 <->/ 2 , 2 , 2\ tr\ 
+L0 l UJ 2 UJ 3 + Ul 2 UJ 3 UO l + UJ 3 00^2 — 2{UJ l + oo 2 + oo 3 ), (5) 

equal to the Ricci tensor component R® which must vanish in the absence of matter 
(vacuum). Nevertheless, we will also discuss cases with 1^0 which represents the 
case of noninteracting matter. 

The question to be settled is the generic behaviour of the system: is it chaotic 
or not? A necessary condition for chaos is the existence of at least one positive Lya- 
punov exponent, but numerical simulations in this direction ^T], [TT], ^| f27], [J7| are 
very difficult, a typical result for the largest exponent being a slightly positive value 
not distinguishable from zero in the vacuum case 1 = 0. More recent computations 



of the correlation dimension [18[ indicate, without definite proof, a probable chaotic 
behaviour. 

The approach adopted here is analytic. The main relevent feature for a global 
knowledge of the behaviour of the system is the singularity structure of the solutions 
[j35[ , p6| . This study, independent of the choice of the metric, must be done in the 
complex plane of r and the issue is: has the system the Painleve property (PP) or 
not? This property is defined as the absence of movable critical singularities in the 
general solution of ([!]), where a singularity is said to be movable if its location in 
the complex plane of r depends on the initial conditions, and critical if branching 
takes place around it (for a tutorial introduction to these questions, see 0). 

The violation of the PP is not enough to decide about chaos. It is in addition 
necessary, but not sufficient, that the general solution takes an infinity of values 
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around the movable critical points |35j (which happens for instance with Log(r— To)). 
On the other hand, the presence of an infinity of movable logarithms is generally 
admitted to be sufficient to ensure the existence of some chaotic regimes. 

Throughout the paper, the term "integrable" means "which has the PP" . 

In section 2, we recall the exact solutions known in closed form. Section |3| proves 
the nonintegrability with a perturbation a la Poincare. In sections 4 and 5, the 
Painleve test isolates all possible particular solutions which can be written as single 
valued expressions in closed form: this proves the inexistence of any vacuum solution 
other than the three known ones. 



2 The three known exact solutions 

The first integral / is identically zero under the condition that the curvature be 



self-dual p] 



a(LogA)' = B + C- A- 2X 1 VBC, X 1 = A 2 A 3 , and cyclically, (6) 
ouj\ = — uji(X 2 uj2 + A3W3), and cyclically, (7) 

where (Ai, A 2 , A 3 ) are constants. The equations ©-© are real only in the Euclidean 
case. The equations for (Ai, A 2 , A 3 ) have three solutions (0, 0, 0), (1, 1, 1), (1, —1, — 1), 
but the third one is not distinct from the second one due to the invariance 
(Ai, A 2 , A 3 , lji, lj 2 , uj 3 ) — > (Ai, — A2, — A3, ujx, — lu 2 , — W3). The two remaining solutions 
define respectively the Euler system (1750) [f| describing the motion of a rigid body 
around its center of gravity 

au[ = and cyclically, (8) 



and the Darboux system ( ||17| eq. (124) p. 149) describing a problem of geometry of 
second degree surfaces 

olJ x = — U1U2 — UJ1UJ3, and cyclically. (9) 

The Euler system has been integrated by Abel and Jacobi (and in the Bianchi 
IX model by Belinskii et ai 0) with elliptic functions, 

A = a(Logpq(A(r - r ), k))', 

B = a(Logrq(\(T-T ),k))', (10) 
C = a(Logxq(A(r - r ), k))', 

depending on the three arbitrary parameters (tq,X, k), where (q,p,r,x) is any per- 
mutation of the letters (s,c,d,n) used in the notation for the Jacobi functions (the 
choice of ref. H] is q=s). The Darboux system has been integrated by Halphen 
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p5| tome I chap. IX) and Bureau [|12j with Hermite modular elliptic func- 
tions. The general solution of the Euler system has for only singularities movable 
simple poles. The general solution of the Darboux-Halphen system is only defined 
inside or outside a movable circle, it is holomorphic in its domain of definition and its 
only singularity is a movable natural boundary defined by the circle. Therefore both 
three-dimensional systems have the PP. These two systems also occur as reductions 



of the self-dual Yang-Mills equations [13, 38 



A third exact solution was found more than forty years ago by Taub [39] who 



noticed that the constraint B 
from six to four 



.4 



C, tu 2 = w 3 is a consistent reduction of the order 



UJ 2 = Co> 3 

B = C = u! 

I 



cosh k\(r 
k 2 



-a 



cosh k 2 {r - r 2 ) 
k\ cosh k\ (r 



n, 



k\ cosh 2 k 2 {r — r 2 ) 



(the above real writing is adapted to a Minkowskian metric a 2 = — 1). It depends on 
the four arbitrary parameters (ti, t 2 , ki, k 2 ), which reduce to three if one satisfies the 
vacuum constraint that I be zero. Its only singularities are movable poles, located 
either at t[ = t% + i(n + 1/2)% /k\ or at t' 2 = t 2 + i(n + l/2)ir/k 2 , with n integer, 
and the four- dimensional dynamical system also has the PP. 



3 Proof of nonintegrability 

The existence of the above three subsystems of order less than six tells nothing 
about the presence or absence of movable branching in the general solution of Bianchi 
IX. 

Following the method initiated by Poincare for the three-body problem 3f|, let 



us look for solutions close to the one (A , B = C ) of Taub in the vacuum case 
a 2 = — 1, k\ = 2k 2 = 2k {k is only defined by its square), in order to introduce the 
three missing arbitrary parameters. One sets 

A = A Q {1 + eAx), B = B (1 + eB,), C = C„(l + eC x ), (12) 

to obtain, at first order in the perturbation variable e, the linear differential system 

A" - 2o~ 2 A\A x = 0, (13) 
Pi - 2a- 2 A B P 1 = Aa- 2 (A B - A 2 )A U (14) 
M'l + 2a- 2 (A B - 2Bl)M 1 = 0, (15) 
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with the notation Pi = Bi + C\, Mi = B\ — G\. The homogeneous parts of ( p~3D and 
( |14[) only differ by the permutation of {k\,T\) and (k2,r 2 ), and the equation (|13|) 
integrates as 

A\ = A; 3 tanh2A;(r — 7~i) + /C4 ^rtanh2A;(r — rx) — — ^ , (£3, k^) arbitrary, (16) 

which only changes the order of the poles of A. The equation fllBD possesses singular 
irregular points of rank two since the coefficient Bq has fourth order poles located 
at the points t 2 . Therefore |28|] its solution contains essential singular points, 
and a sufficient condition that they be critical for ([15]) is that either the rank of the 
singular irregular points be odd (which is not the case here), or the formal solution 
in the neighborhood of t 2 

+00 +00 
Mi = k 5 e a +^-rD £ A + (r _ T tf+s + + he «-/(r-r D J2 A-( T - r' 2 ) k+s -, A + A ^ 0, 

k=0 k=0 

(17) 

where (a±, s±, A fc ) are two sets of constants, have one of its two Thome indices s± 
irrational. The two values (a±,s±) are given by the first two most singular terms 

a ± = ±k^ smh2k(r 2 - r[), s± = 1 =f 2 cosh 2k(r' 2 - t[). (18) 

The generically complex values, with irrational real and imaginary parts, for the 
Thome indices s± introduce a nonalgebraic (i.e. transcendental) critical branching in 
Mi and therefore a movable critical transcendental essential singularity in the general 
solution of the Bianchi IX model. It may be useful to insist on the impossibility to 
remove this critical singularity by considering either Stokes sectors or some algebraic 
transform of (A, B, C), as would be the case with u 2 for the equation whose general 
solution is 



u = Kie T ~ T o ■sjr — r + K 2 e T ~ T o jy/r — r , (K x , K 2 ) arbitrary. (19) 

The reason for this impossibility is the irrational nature of the two Thome indices. 

Thus, the Bianchi IX model has not the Painleve property and its general solu- 
tion contains an infinity of movable logarithms. As explained in the introduction, 
it is therefore quite probably chaotic, a fact supported by some recent numerical 
simulations 



In the limit ki — > 0, k% — ► where the solution of Taub degenerates into the 
two-parameter solution ([TJJ p. 150) of the Darboux-Halphen system 

A 1 B C t — r\ .. it 

- = 7- — = — = 7Z ITT?' k n T ± ~ r i) = fc 2(T2 - t 2 ) = 1-, (20) 

a r — r[ o o [t — T 2 y I 

the movable essential singularity becomes noncritical (a± = ^2{t 2 —t[),s± = 1 =F 2) 
and the formal series for Mi can be summed with the Euler function E(x) 

Mi = k^x' 1 + k 6 e~^(l - x + 2x 2 - E(x)), x = J ~ ^ (21) 
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E{X)= I- —u d "- (22 > 



It could be interesting to consider an equation equivalent to (15) with polynomial 



coefficients, by the change of variables T = tanh c(r — r 2 ), M\ = mi cosh c(r — r 2 ): 
4a 2 (l-T 2 ) 2 ^ 

+ [a 4 (l + Tf + (1 - T) 4 + 2a 2 (l - T 2 ) 2 + 8a 2 T 2 - 4a 2 ]m! = 0, 

a = e 2c(T2 " Tl) , (23) 

which sends the singular irregular point to T = oo. 

Remark. A perturbation of the Taub solution has already been performed || up 
to and including second order, without finding critical singularities. The effect is 
due to the neglection of terms involving A in the three right hand sides of (|l|) (see 
system (4.4)-(4.6) in ref. M), and the present results a posteriori prove that these 
terms are crucial. 



4 Local singularity analysis 

Let us now try to obtain the list of all the exact solutions which can be written 
in closed form. Such solutions can only be particular and depend on at most four 
arbitrary parameters, since they cannot |35| contain any of the two infinities of 



movable logarithms found in section |3|. Three of them are already known, see section 
2. To achieve this goal, let us express necessary conditions for the absence of any 
infinity of movable logarithms in local representations of the general solution. The 
local, necessary information thus obtained will then have to be transformed into a 
global, sufficient one in the shape of a closed form solution. The method involved is 
called singularity analysis ||. 

The singularity analysis of this dynamical system has recently been undertaken 



[ Hp and there exist two different local representations of (A, B,C) by Laurent series 
bounded from below, i.e. describing a meromorphic-like behaviour. 

The first family of movable singularities, about a point which we denote r±, is 

A/a = x' 1 + a 2 x + 0(x 3 ), X = r - n, 

B/a = b x + b lX 2 + O( X ;i ), (24) 

C/a = c x + c lX 2 + O(x 3 ), 

I /a 2 = 6a 2 - 2(b + c ) + &iCi/(& c ), 

and its Fuchs indices are (—1, 0, 0, 1, 1, 2), corresponding to the orders at which the 
six arbitrary coefficients enter the expansion, respectively (r 1; b , c , &i, Ci, a 2 ). 
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The second family, about a movable singularity which we denote r 2 , is 

A/ a = x' 1 + «2X + 0( X 3 ), X = t-t 2 , 

B/a = x~ 1 + b 2 x + 0(x 3 ), (25) 

C/a = X - l + c 2 x + 0(x 3 ), 

I jo 2 = -6(a 2 + b 2 + c 2 ), 

its Fuchs indices are (—1, —1, —1, 2, 2, 2), but the arbitrary coefficients (r 2 , a 2 , & 2 , c 2 ) 
correspond to only four Fuchs indices, respectively (—1,2,2,2). 

The first series represents a locally meromorphic behaviour of the general solu- 
tion. On the contrary, the second series only represents a four-parameter locally 
single valued particular solution, and it tells nothing about some possible multival- 
uedness about r 2 in the general solution. The reason is the presence of two negative 
indices in addition to the ever present —1, counted once even if it is multiple (con- 
trary to the erroneous argument presented in |15|). In such a situation, the method 
of pole-like expansions initiated by Sonia Kowalevskaya p0| , |31| , 22 , [L] is unable to 



represent the general solution and to check the absence of movable branching about 

72- 

Therefore, one must perturb the second Laurent series in order to extract the 
information contained in the two remaining negative indices. 

The similar local representations of the particular solutions are the following. 
For the Belinskii et al. solution ([II]) of the subsystem (|8p, the family ( |2"1| ) with bi = 
Ci = 0, 3a 2 = b + c and indices (—1, 0, 0), and the family (|25|) with a 2 + b 2 + c 2 = 
and indices (—1,2,2). For the Halphen solution of the Darboux subsystem (||), 
the family fl2j|) with bi = 2& y / c ~, c\ = 2c a/5o", a 2 = (\/&o — v /c ") 2 /3 and indices 
(— 1, 0, 0), and the family 

A/a = x' 1 + a-iX -2 + • • • , X = t ~ t 2 , 

B/a = X' 1 + b_ lX - 2 + ..., (26) 

C/a = x' 1 + c-ix~ 2 + • • • , 

with indices ( — 1,-1,-1). This "descending" Laurent series shares its first term 
with ([25]) and it evidently describes a solution missing in (|25|). 

As to the Taub particular solution of the full system (|l|), near t[ it is represented 
by (^|) with b = c ,bi = c x but, near t' 2 , it is impossible to fit its Laurent series 
into one of the above ones. This is due to the presence of the movable essential 
singularity detected in perturbation, which makes meaningless the computation of 
Fuchs indices (|28| chap. XVII). Its Laurent series near t' 2 will be obtained in next 
section. 

5 All possible solutions in closed form 

An algorithmic method to perform the perturbation of the second family (|25| ) has 
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recently been proposed [^(], [H] , and later shown to be a natural application of the 
theorem of perturbations of Poincare and Lyapunov. At each perturbative order, 
this method checks the absence of movable logarithms likely to occur at all indices, 
then adds to the previous Laurent series another Laurent series, also bounded from 
below but starting with a strictly lower singularity power. The resulting "doubly 
infinite" Laurent series, which is the well known standard local representation of 



a single valued function, is the union of the two pieces (|25f) and (p6|) , with some 
"interference terms" of the highest interest. 

Our interest here is not to prove multivaluedness, which has been done in section 
|3|, but to exploit the necessary conditions for the absence of movable logarithms in 
order to detect particular locally single valued solutions, which will then have to be 
written in closed form if the conditions are sufficient. 

In our case where the smallest index is —1 and the highest one 2, to perform this 
perturbative test up to some maximal order N, one sets \TA 

A N 2+N-n 

— = X~ 1 J2 £n E af x\x = T--T 2 , and cyclically, (27) 

a n=0 j=-n 

and puts these polynomials (with negative and positive powers of x) hito a polyno- 
mial definition of the dynamical system 

a A E A = a 2 (AA" - A' 2 ) - A* + A 2 {B - C) 2 = and cyclically. (28) 

Writing the LHS E A of fl28|) as a polynomial similar to (|2T|) 



N 2+N-n 

E A = X~ 4 J2 £n E E A 3 X 3 + ••• and cyclically, (29) 

n=0 j=—n 

one then solves the set of equations = by respecting the ordering that E^? ■* = 

must be solved before E^ = 0, with n! < n,j' + n' < j + n, (n',j') ^ (n,j). The 
first equation (n,j) = (0,0) is the only nonlinear one 

- 4 0)2 + (b { o ] - 4 0) ) 2 1 = and cyclically, (30) 



F - 



and it has already been solved, see eq. ( p5|) or (p6]), as = = Cq^ = 1 since a is 
only defined by its square. As usual for a perturbative method, all other equations 
are linear 



V(n,j)^(0,0): P(j) 



b {n) 



Qf=0, P( J ) = -( J + l)( J -2)J, (31) 



where P is a multiple of the identity matrix X independent of n and the column 
if 



vector only depends on previously computed coefficients. 



S 



At every perturbation order n and every triple Fuchs index j G {—1,2}, the 
linear system for , fej , is singular and has a zero rank; if the orthogonality 

condition = is satisfied, three more arbitrary coefficients enter the expansion, 
but if it does not an impossibility occurs and logarithms must be introduced. In 
the former case, the arbitrary coefficients need only be introduced the first time 
the Fuchs index is encountered, i.e. (n,j) = (1, —1), (0,2), because next times they 
would only perturbate these ones. Since the formal solution evidently does not 
depend on the seven parameters (r 2 , , £0,-1, sb_\, scl{), one can freely 

choose one of the values of the four quantities (r 2 , cll{ } b_\, cl{) related to the triple 
Fuchs index —1. This choice of gauge will be used below. 

Due to the nice symmetry of the second family, only one of the three scalar 
equations of (|3~lf) needs to be considered, the results for the two others being derived 
by cyclic permutations. Parity makes hence aj 1 vanish for n + j odd. 

The costless violation occurs at third order for the Fuchs index —1 



(0) 


= 1, a_\ = arbitrary, = arbitrary, 


(32) 


4a« 


= 4aL?-(^-a 2 , 


(33) 


4a® 


= 4ag 3 -(ag + a + c3)(6g-c«) 2 , 


(34) 


4aS° 


= (6 W _ c «){(54 0) - 46f - 44 0) )(foL 1 l - c ( H) 






-(2aL 1 l-a-cS)(4° ) -4° ) )} 


(35) 


a? 


= -aPaM + (bP-c?Wi-£l), 


(36) 




= (3, — 1) : violation. 


(37) 



Unless the three conditions Q 



(3) 



E 



(3) 



are satisfied, movable logarithms 



start entering the expansion and the Painleve test fails. The formal solution will 
thus contain an infinity of logarithmic terms. The three conditions Q 



(3) 



are 



not independent and, together with the first integral /, they are best expressed in 
coordinates adapted to the ternary symmetry: 





R 



i 



9Ci, Q£ 



(3) 
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-i 



2 2 
C\ = —x{ux 2 — Avxy + uy 2 ) = 

C 2 = vx 3 — 2ux 2 y + vxy 2 + 2uy 3 = 0, 

2 - -6y/3w + (18 terms)e 2 + 0(e 4 ) 



R 



Qt 



o 

-V2/V3 
1/V3 
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1/V2 
1/V6 
1/V3 

9V3„ 

- L-2 



-1/V2\ 

i/VQ 

VV3 y 



I/a 
Equations C\ 



(38) 

(39) 
(40) 
(41) 



0, C*2 = are independent of (w, z) and have five solutions 

: bf = 4 0) , 6« = c« (42) 



u 



x 
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u = v = 0: 4 0) = 4 0) = 4 0) , (43) 
x = y = : aL 1 } = = c% (44) 

corresponding to only three distinct cases. 

The first constraint (|42|) implies the equality of two of the components (A, B, C) 
at every order and thus represents the four-parameter solution of Taub ([11]) . 

The second constraint ([f3|) is not refined at next condition (n,j) = (2,2) but it 
is restricted at (n,j) = (5, —1) by the relation wx(x 2 — 3y 2 ) = 0, which splits into 
either x[x 2 — 3y 2 ) = 0, i.e. a subset of case (|42|), or w = 0, so the second constraint 



( |43"D becomes u = v = w = 0. This represents the three-parameter solution of the 
Darboux-Halphen system (Q). 

For the third and last constraint (fill), the doubly infinite Laurent series (p7|) has 
the same sum than the semi-infinite Laurent series Q2"5D because one can assign the 
gauge a_\ = = c_] to zero, exactly like the series of negative powers J2jl-oo( T 2 ~ 
To)~- , ' -1 (t — r o y , convergent outside a disk centered at r containing r 2 , and the series 
of positive powers — I]j=o( r 2 — r i)' ?_1 (' r — r i) : '5 convergent inside a disk centered at 
Ti not containing T2, both represent the same function (r — T2) _1 . The perturbative 
Painleve test reduces in this case to the test of Kowalevskaya and one cannot rule out 
the possibility of a four-parameter, global, closed form, single valued exact solution 



extending the three-parameter solution (|iOD of Belinskii et al, which it contains for 



02 + &2 + c< 2^ — V^w = 0. We have not yet succeeded in finding this closed form. 

Since the local singularity analysis shows that the only movable branching is 
logarithmic (no algebraic branching such as x q i ( l rational), this solves entirely the 
question of closed form global solutions in the vacuum case 1 = 0: there exists no 
solution to the Bianchi IX model in vacuum other than the three known ones. Let 
us recall once more that the local, not in closed form solutions are illusory [[J^j . 

Remarks. 

1. The lowest perturbation order (two) at which the failure of the Painleve test 
occurs is one more than the perturbation order of the solution of Taub, be- 
cause the present perturbation starts from a formal solution with a simple, 
not double, pole. 

2. This dynamical system is one more example |14[ of a family of movable sin- 



gularities with negative Fuchs indices which contains all the information on 
the integrability, while the family with positive indices, in this particular case, 
contains no useful information at all. 
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6 Conclusion 



With two different proofs using perturbative methods, we have shown the non- 
integrability in the Painleve sense of the Bianchi IX model, by exhibiting mov- 
able essential transcendental critical singularities which quite probably imply chaos. 
Moreover, we have proven the inexistence of any vacuum solution other than the 
three known ones. 

One important open problem is the understanding of the apparently zero value 
for the highest Lyapunov exponent. Has this something to do with the existence of a 
local singlevalued formal representation of the general solution by the meromorphic 
series (21), whose radius of convergence could then be worth studied? 

Our result should probably not affect the overall understanding of the Kasner 
epochs, for which a detailed description has been given || |32|, |?], [l9| . 
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Note added in proof. After the submission of this paper, we obtained a copy of a 
paper |16|] correcting a previous one |L5 1 : the authors use the perturbative Painleve 



test for testing the negative indices and find movable logarithms. 
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